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ABSTRACT 

We present a model of a dusty disk with an inner hole which accounts for 
the Spitzer Space Telescope Infrared Spectrograph observations of the low-mass 
pre-main sequence star CoKu Tau/4. We have modeled the mid-IR spectrum 
(between 8 and 25 /xm) as arising from the inner wall of a disk. Our model disk 
has an evacuated inner zone of radius ~ 10 AU, with a dusty inner "wall", of 
half-height ~ 2 AU, that is illuminated at normal incidence by the central star. 
The radiative equilibrium temperature decreases from the inner disk edge out- 
ward through the optically-thick disk; this temperature gradient is responsible 
for the emission of the silicate bands at 10 and 20 fim. The observed spectrum 
is consistent with being produced by Fe-Mg amorphous glassy olivine and/or py- 
roxene, with no evidence of a crystalline component. The mid-infrared spectrum 
of CoKu Tau/4 is reminiscent of that of the much older star TW Hya, where it 
has been suggested that the significant clearing of its inner disk is due to planet 
formation. However, no inner disk remains in CoKu Tau/4, consistent with the 
star being a weak-emission (non-accreting) T Tauri star. The relative youth of 
CoKu Tau/4 (~ 1 Myr) may indicate much more rapid planet formation than 
typically assumed. 
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Subject headings: planetary systems: formation, protoplanetary disks, stars: pre- 
main sequence, circumstellar matter 

1. Introduction 

The formation of planets around young stars is thought to involve the eventual clearing of 
the natal dusty protoplanetary disks, as planets accrete and sweep out material. Coagulation 
and accretion processes are expected to occur fastest in inner disk regions, where surface 
densities arc higher and orbital periods shorter (e.g., Wuchterl, Guillot, & Lissauer 2000, 
and references therein). An observational signature of inner disk gaps or holes would be 
the reduction or elimination of emission from hot dust, corresponding to a decrease in the 
short-wavelength infrared emission. 

The low-mass pre- main sequence star CoKu Tau/4 exhibits a spectral energy distri- 
bution (SED) with properties suggesting a substantial dusty disk with an evacuated inner 
region; strong infrared excess emission is detected at wavelengths > 10 fim, with little or no 
excess over photospheric fluxes at shorter wavelengths, as shown by observations with the 
Infrared Spectrograph (IRS) on the Spitzer Space Telescope (Forrest et al. 2004). The inter- 
pretation is that CoKu Tau/4 has a dusty disk which is mostly cleared of small dust grains 
inside a radius of 10 AU. The quahtative nature of the SED of CoKu Tau/4 is reminiscent 
of that of the Classical T Tauri stars TW Hya (Calvet et al. 2002; Uchida et al. 2004), CM 
Aur (Marsh & Mahoney 1992, 1993; Rice et al. 2003), and the Herbig Ae star HD 100546 
(Bouwman et al. 2003). Quillen et al. (2004) suggest that the inner disk hole in CoKu Tau/4 
might be produced by a relatively low- mass planet; if this is the case, CoKu Tau/4 is par- 
ticularly interesting because its relative youth (~ 1 Myr vs. 10 Myr for TW Hya; Webb et 
al. 1999), places severe constraints on planetary formation (e.g., Wuchterl et al. 2000, and 
references therein). 

In this paper we present radiative transfer models of CoKu Tau/4 which account for 
the IRS spectrum and other constraints on the SED. We confirm our previous estimate 
of the disk inner radius (Forrest et al. 2004), and show that the observed 10 /xm silicate 
emission feature is likely to be produced in the hot outer layer of the disk "wall" facing the 
star, without requiring any small dust inside the inner hole. We also comment on the dust 
compositions consistent with the observations. 
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2. Stellair properties, extinction, and SED 

CoKu Tau/4 is a weak-emission T Tauri star (WTTS) of spectral type of Ml. 5 (Cohen 
& Kuhi 1979), with no detectable short wavelength (UV) excess and an Ha equivalent width 
~ 1.8 A (Cohen & Kuhi 1979). We assume a distance of 140 pc to the star, based on its 
location in the Taurus- Auriga molecular cloud (Kenyon, Dobrzycka, & Hartmann 1994). 

To estimate the stellar parameters, we need to quantify the extinction to the star. Cohen 
& Kuhi (1979) found Ay = 0.94 ± 0.29 from their spectral type determination, adoption 
of the standard reddening law, and using fluxes from their scanner spectra in wavelength 
regions near 5400, 6000, and 6700 A. However, such a low extinction is not consistent with 
relatively recent observations of near-infrared colors. Kenyon & Hartmann (1995; KH95) 
flnd J - H = 1.10, and the 2MASS catalog yields J - H = 1.08, while from Table A4 of 
KH95, the unreddened color of an Ml. 5 star should be J — H = 0.67. Using an extinction 
law with Ry ~ 3.6, as indicated by the recent study of Calvet et al. (2004), this would 
suggest Av ~ 3. The H — K colors suggest slightly larger extinctions, but the latter result 
is much more sensitive to errors in photometry and extinction laws than the result obtained 
from the J — H colors. 

Recent optical spectra of CoKu Tau/4 from the FAST spectrograph on the Fred L. 
Whipple 1.5m telescope on Mt. Hopkins (discussed in Kenyon et al. 1998) were inspected. 
Although conditions were not photometric, comparison with other stars observed on the 
same night indicate a much redder spectrum than shown by Cohen and Kuhi (1979), roughly 
consistent with Ay ~ 3, in agreement with the values of extinction inferred from 2MASS 
colors. For the purposes of this paper we therefore adopt Ay — 3. Further optical monitoring 
of CoKu Tau/4 would be desirable. 

Following KH95, we assign a stellar effective temperature of T* =3720 K to CoKu Tau/4. 
For Ay = 3 we find a stellar luminosity of approximately = 0.61 Lq, with an error of 
order 20% or more given possible variability and the uncertainty in extinction corrections. 
(According to 2MASS, J = 10.16, while KH95 cite J = 10.37.) Evolutionary tracks from 
Siess, Dufour, & Forestini (2000) give a mass of ~ 0.5 Mq and the age of ~ 1.2 Myrs for 
the central star, while those of Baraffe et al. (1998), give 0.66 M0 and 1.6 Myrs, respectively. 
Both ages are consistent with that of the rest of the population of Taurus, ~ 1 - 2 Myrs 
(Hartmann 2003). 

There are two alternatives for the cause of the high extinction observed towards CoKu Tau/4, 
either of interstellar/intercloud origin, or, if CoKu Tau/4 were an edge-on disk, of local ori- 
gin. However, edge-on disks are underluminous, because the direct light from the star gets 
heavily extinguished, so the optical spectrum is stellar scattered light (D'Alessio et al. 1999). 
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This scattered component is fainter than the stellar direct flux because the scattering re- 
gion does not cover 4 tt steradians. A correction of Ay = 3 applied to this scattered light 
component would still leave the star under luminous; however, with Ay = 3 and its current 
IR magnitudes, the properties of CoKu Tau/4 agree with those of the low-mass population 
of Taurus. In addition, the large constrast between 10 and 20 fim in CoKu Tau/4 is much 
larger than predicted in high inclination disks with emission in 10 fim (Pascucci et al. 2004) 

We then adopt the hypothesis that the extinction towards CoKu Tau/4 is interstel- 
lar /intercloud; the presence of a scattered light nebulosity ~ 1' from CoKu Tau/4 (D. Pad- 
gett, unpublished) supports the assumption that it is located in a high density region of 
the molecular cloud. With this assumption, the IRS spectrum must also be corrected for 
extinction. Figure 1 shows the IRS spectrum corrected using two different reddening laws, 
one proposed by Draine (2003) and the other one by Moneti et al. (2001) based on ob- 
servations from Lutz (1999) and 10 and 18 //m profiles from Simpson (1991). Between 1.2 
and 8 /xm, the reddening law proposed by Moneti et al. (2001) agrees, within the uncer- 
tainties, with the reddening law recently determined using IRAC data (Indcbetouw et al. 
2004). The overall features of the spectrum shown in Forrest et al. (2004) are preserved, 
namely, a steep, roughly photospheric spectrum shortward of about 8 fim, the weak silicate 
feature at 10 /im, the steep rise between 14 and 20 /xm, and the flattening of the emission 
to longer wavelengths. The shape of the 10 //m sihcate feature is slightly dependent upon 
the extinction law used. The flux error at each wavelength was derived by taking half of the 
absolute value of the difference in flux from the spectra independently obtained at the two 
nod positions of the telescope for each order of Short- Low (SL; 5.2-14 /xm, A/AA ~ 90) and 
Long- Low (LL; 14-38 //m, A/AA ~ 90) (See Forrest et al. 2004). 

Figure 2 also shows fluxes from 2MASS and long-wavelength measurements from IRAS 
from Strom et al. (1989), which may be contaminated by extended emission. For comparison 
we show a synthetic stellar spectrum from Bruzual & Chariot (1993) for T*=3720 K, — 

1.9 R0, and = 0.5 M0, along with the median SED of Classical T Tauri stars in the 
Taurus Association (D'Alessio et al. 1999), scaled to the flux in the H band. It is evident 
that the mid infrared emission of CoKu Tau 4 is significantly less than those of typical 
Taurus stars whose disks lack inner holes. In contrast, the flux levels in the 20-35 /^m region 
are relatively bright compared with the median Taurus SED. Our task is to explain this 
spectrum. 
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3. Model 

We have constructed models for CoKu Tau/4 that are similar to the one calculated for 
TW Hya (Calvet et al. 2002; Uchida et al. 2004). The basic model for TW Hya is that of a 
fairly massive Classical T Tauri disk truncated at an inner radius of about 4 AU (Calvet et 
al. 2002); at this radius the optically-thick disk is frontally-illuminated by the central star. 
A major difference between TW Hya and CoKu Tau/4 is that the former object exhibits 
clear evidence for gas accretion onto the star (MuzeroUe et al. 2000; Batalha et al. 2002). 
The SED of TW Hya also exhibits a substantial excess at 10 fim and even some small excess 
emission at wavelengths as short as ~ 5 /im (Sitko, Lynch, & Russell 2000). Calvet et al. 
(2002) interpreted these results to mean that there is a small amount of optically-thin dust 
and gas interior to the optically-thick disk wall, which accounts for the hotter dust emission 
and the accreting gas. In contrast, CoKu Tau/4 is a WTTS, with a small H^ equivalent 
width consistent with stellar chromospheric emission; the accretion rate in this system must 
be extremely small, if present at all (MuzeroUe et al. 2000). In addition, the SED of 
CoKu Tau/4 is consistent with purely stellar photospheric emission at wavelengths A < 8 
//m (Forrest et al. 2004). As discussed in the previous section, there are uncertainties at the 
~ 20% level in extrapolating the photospheric continuum from shorter wavelengths to the 
IRS range; however, the spectral slope at the short-wavelength end of the IRS spectrum is 
clearly consistent with being purely photospheric, as are IRAC colors of this object between 
3.6 and 8 /xm (Hartmann 2004, personal communication). 

Thus, we model CoKu Tau/4 as a central star with the properties described in §2 
surrounded by an optically-thick, truncated disk, with an inner wall at a radius -Rwaii- Inside 
this radius the disk is optically thin; in §3.4 we estimate the maximum amount of dust in the 
inner hole or gap consistent with the spectrum. Outside i?waii there is a disk which might 
contribute to the SED at wavelengths longer than 25 //m. At R — -Rwaii, the disk receives 
radiation frontally from the central star; as we detail in the following discussion, the wall 
appears to be responsible for most, if not all, of the excess emission detected by IRS. 

3.1. Treatment of the inner disk "wall" 

The wall is the transition between the inner optically thin disk or gap and the outer 
disk. We assume that it is uniform in the vertical direction and solve its radial structure. 
Closer to the star, the wall has an optically thin atmosphere, and its radial optical depth 
increases with radius. Figure 3 shows schematically the structure of this wall. 

The temperature distribution of the wall is calculated using the procedure described by 
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Calvet et al. (1991, hereafter CPMD; 1992), which has been successfully apphed to calculate 
the radial structure of the wall at the dust sublimation radius in disks around Classical T 
Tauri Stars by MuzeroUe et al. (2003). Since we relax some approximations made in previous 
papers, we write the equations in detail. 

We assume that the stellar radiation penetrates the wall with an angle zero between the 
incident direction and the normal to the wall's surface, and that the radial thickness of its 
atmosphere is negligible compared to its distance to the central star. The incident radiation 
field is separated into two wavelength ranges, one characteristic of the disk (related quan- 
tities have subindex "d") and one characteristic of the incident stellar radiation (quantities 
with subindex "s"). In order to find an analytical solution for the wall radial distribution 
of temperature, we assume that the opacities are constant in the wall atmosphere and that 
scattering of stellar radiation is isotropic. We also assume that there is no heating source 
other than the incident and scattered stellar radiation. At every depth into the wall at- 
mosphere the net outward radial flux at the disk frequency range is equal to the absorbed 
radiative flux at the stellar frequency range (see CPMD), i.e, the Eddington flux , Hd, is 
given by 

Hairs) = ^a[{l + Oe'^' + C^e'^'^] (1) 

where a = 1 — w, w is the mean albedo to the stellar radiation [w = (Js/Xs)-, P = (Sa)^/^, 
Fo = L*/47ri?^jjii and Tg is the total mean optical depth to the stellar radiation in the incident 
direction (i.e., along rays parallel to the disk mid- plane, perpendicular to a cylindrical surface 
at R — i?waii)- The constants in equation (1) are given by 

' /3[1 + (2/3/3)] (1-/32)- ^' 
The zeroth moment of the transfer equation can be written as (Mihalas 1978), 

f^-J,-Sa (4) 

where is the total mean optical depth at the disk frequency range, Sd is the local source 
function and Jd is the mean intensity. CPMD assume strict LTE, i.e., Sd = Bd = ctrT^/tt, 
where is the Stefan-Boltzmann constant. Here, we include a term corresponding to the 
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emissivity by scattering at the disk frequency range, i.e., 

Sd = {ndBd + (^dJd)/{f<'d + o-d), (5) 

assuming isotropic scattering (Mihalas 1978). With this source function, the zeroth moment 
of the transfer equation is given by 

Bd^Jd (6) 

i^d dTd 

where Xd = i^d + CTd- 

Using the Eddington approximation for the disk radiation field, the first moment of the 
transfer equation is 

3//, (7) 



dTd 

where we have defined q = Xs/Xd- Integrating Jd from equation (7), using Jd(0) = 2Hd{0) 
as the boundary condition, and substituting it along with dHd/dTd from equation (1) into 
equation (6), the temperature as a function of can be written as 



where 



Td{Tdf = a^[C[ + C'.e-'^^- + C'^e-^'^^^] (8) 



C; = (l + C0(2+?)+C.(2 + ±) (9) 



q) \ (3q. 

C^ = (l + Ci)f^--) (10) 



The temperature at the innermost radius of the wall, at the top of the wall atmosphere 
where Td — Tg — 0, is Tq, can be used to define the wall's radius. 



-Rwall — 



(1 + C) (2 + q^) +CJ2 + Pq^)] (12) 



1/2 I 



We use To as an input parameter to characterize the atmosphere of the wall. The value of 
-Rwaii follows from Tq and the adopted dust properties, using equation (12). The mean opacity 
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in the disk wavelength range is calculated using the Planck function evaluated at Tq as the 
weighting function, and the mean opacity in the stellar wavelength range is calculated using 
the Planck function at T^. We use Planck means in order to have the correct temperature 
in the optically thin layers of the wall atmosphere (Mihalas 1978). 

The minimum temperature of the wall is obtained when — > oo and it is given by 



T 



mm 



Figure 4 shows the temperature as a function of the disk mean optical depth (measured 
from the innermost radius of the wall toward larger radii) for two models with the same value 
of Tq but different kind of sihcate: Mgo.s Feo.5 olivine and Mgo.s Feo.2 pyroxene. Both models, 
with appropriated values for the disk inclination angle and the wall height, reproduce the 
observed spectrum of CoKu Tau/4(see §3.3). The temperature of the wall as a function 
of optical depth depends on the monochromatic absorption coefficient of the dust through 
the quantity q, so it has to be calculated for every dust composition attempted to fit the 
spectrum. However, for a given composition, the temperature is independent of the total 
dust to gas mass ratio, since this ratio affects the heating and the cooling of the grains in 
the same way, canceling out. 

We have calculated the asymmetry factor g for Mg Fe olivine and pyroxene, the materials 
that provide the best fit to the observed spectrum (sec §3.3), assuming grains are compact 
spheres and using the method developed by Wiscombe (1979). We note that the asymmetry 
factor g of these ingredients decrease with wavelength. In particular, for a maximum grain 
size Qmax — 0.25 /xm, g ~ 0.6 for A = 0.55 //m and g ~ 0.2 for A = 1.22 //m. Taking a Planck 
mean of the g factor using the Planck function evaluated at the stellar effective temperature 
as the weighting function, the mean asymmetry factor for scattering of stellar radiation is 
< g >~ 0.3. This justifies the assumption of isotropic scattering {g ~ 0) we have used to 
find an analytical solution for the temperature distribution in the disk. 



3.2. Spectrum of the wall 

Accounting for the limb brightening produced by the temperature inversion in the wall 
atmosphere, we calculate the emergent intensity for each surface element of the wall as seen 
in the plane of the sky. The geometry is described in detail in the Appendix. Here we briefly 
outline the required expressions. 

The cosine of the angle between the normal to the surface of the wall and the fine of 
sight is given by 
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cos © = sin i sin 9 (14) 

where i is the angle between the disk rotation axis and the hne of sight and 6 is such that the 
coordinate X in the plane of the sky is X = -Rwaii cos 9. Each element in the visible surface 
of the wall has a thermal emergent intensity given by 

1 

L = / B,[Ta{Ta)]e-^^l"''^dT, (15) 

cos k) Jq 

where = t^^k,^ / Xd) ■ The emergent flux is calculated by integrating the emergent intensity 
over solid angle on the visible surface of the wall. The limits for the integral are given in the 
Appendix. 

Using the intensity given by the equation (15) we obtain a limb brightening effect 
because rays that pierce the limb of the wall (i.e., x ~ Rwaiu see Figure 9) have an optical 
depth equal unity in layers with higher temperature than rays that intersect the wall closer 
to the center (i.e., a; ~ 0). This effect would be neglected if the thermal emergent intensity 
is approximated as isotropic and evaluated at cos = 1 with the expression 



OO 



B,[Ta{Ta)]e-^''dT,. (16) 



In this approximation, the emergent flux is calculated by multiplying the total intensity 
by the total solid angle subtended by the visible side of the wall, which at a distance d from 
the observer is given by 



Q^^u = 2 cosi {^^^ (5VT^ + sin-^ 5) if 5<l (17) 



^wall = TT ( 1 cosi if 5 >1 

5^:^tanz, (19) 

and i?waii the (vertical) height of the wall photosphere (measured from the midplane, in the 
vertical direction). Equations (17) and (18) are consistent with the sohd angle portion of 
the expressions used by DuUemond, Dominik & Natta (2001) to calculate the flux of an 
isothermal wall. 



where 
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Limb brightening effects are important when trying to fit a high resolution spectrum 
as those of Spitzer-IRS. For instance, we find that taking into account hmb brightening 
the temperature Tq required to fit the spectrum is smaller by 5 K than the one in the 
approximated solution, increasing the wall radius by 10 %. In turn, this translates into a 
height of the wall 20% smaller than the one inferred for the approximated solution. Even 
more important, the shape of the 20 iJ,m region of the spectrum only can be fitted with the 
detailed model; the approximated model shows an emission deficit for all the compositions we 
have tried (described in §3.3). The 20 /im flux deficit of the approximated model compared 
to the observed spectrum could be wrongly interpreted as the lack of additional unknown 
opacity sources. 



3.3. Dust properties 

We use a mixture of grains composed of silicate with mass fraction (gu =0.0034, organics 
with (org =0.001, and troilite with Qroi — 7.68 x 10~^. The grains are assumed to be 
spheres. We adopt the standard MRN grain size distribution n{a) ~ a"^-^ (Mathis, Rumpl 
& Nordsieck 1977), between minimum radius amin and maximum radius amax- The opacity 
is calculated using the Mie theory with a code developed by Wiscombe (1979). Optical 
constants for the organics are taken from Pollack et al. (1994, hereafter P94) and for troilite 
from P94 and Begcmann ct al. (1994). Sublimation temperatures Tgn = 1400 K, T^rg = 425 
K, and Ttroi = 680 K are adopted. 

The shape of the 10 fim silicate band and the position of the peak in the spectrum 
of CoKu Tau/4 indicates that the grains responsible for the emission are small. We adopt 
o-min = 0.005 fim and a^nax = 0.25 fim, which are consistent with ISM grains (Draine & Lcc 
1984). If grains were larger than this, then the peak of the silicate emission would broaden 
to longer wavelengths than observed. 

We have considered several possibilities for the composition of the silicate (sec Hcnning 
et al. 1999^). The observed 10 fim band is smooth and narrow, with no evidence of the 
substructure characteristic of crystalline components such as enstatite, forsterite and silica 
(e.g., Uchida et al. 2004), suggesting the silicate should be glassy and amorphous. In 
particular, the shape of the observed 10 /im band allows us to rule out glassy bronzite (with 
optical constants from Dorschner et al. 1988), crystaUine bronzite (optical constants from 
Hcnning & Mutschke 1997), crystalline silicate of mean cosmic composition Mgo.50 Feo.43 
Cao.03 Alo.04 Si03 (optical constants from Jager et al. 1994), oxygen rich and oxygen deficient 



http : //www . astro . uni- j ena . de/Laboratory/Database/ jpdoc/index . html 
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circumstellar silicate (optical constants from Ossenkopf et al. 1992), and the absorption 
coefficient proposed by Draine (2003^ ) for interstellar medium dust. The Draine(2003) 
mixture and the oxygen rich and deficient circumstellar silicate from Ossenkopf et al. (1992) 
produce an excess flux at the 10 - 13 fim wavelength range and a deficit at longer wavelengths 
compared to the CoKu Tau/4 spectrum. The model spectrum calculated for glassy silicate 
of mean cosmic composition Mgo.50 Feo.43 Cao.03 AI0.04 SiOa (optical constants from Jager et 
al. 1994) fits well the 10 //m band, but the shape of the spectrum between 20 and 30 //m 
is different from the observed one. Figure 5 shows some examples of the synthetic spectra 
of models that we have considered do not fit the observed spectrum. It is important to 
mention that we have changed the temperature Tq, inclination angle and height of the wall 
in the ranges described below and we have not found any combination of parameters for 
these compositions that produces a good fit to the observed spectrum. 

On the other hand, amorphous Mg-Fe glassy olivine ( Mgo.5 Feo.5 Si04) and glassy 
pyroxene (in particular, Mgo.95 Feo.05 SiOa and Mgo.s Feo.2 SiOa) with optical constants from 
Dorschncr ct al. (1995), produce a reasonable fit to the spectrum. In particular, the pyroxene 
has a broader 10/xm band than the observed one but fits the spectrum for A > 18 /im better 
than the olivine. On the other hand, the shape of the 10 fim band predicted for the olivine 
is in better agreement with the observed spectrum. We find that if olivine or pyroxene has 
less magnesium then the resulting model spectrum has a flatter 10 /im band compared to 
the observed one. However, we are not considering both olivine and pyroxene, nor the effect 
of the grains shape and porosity, which could affect the shape of the bands. Thus, given the 
errors and our limited set of compositions, we only conclude that both, glassy olivine and 
glassy pyroxene can be present in the wall of CoKu Tau/4, and that there is no evidence of 
a crystalline component. Figures 6 and 7 show examples of the model spectra calculated for 
each of these ingredients that flt reasonably well the observed spectrum. 

For each grain composition we calculate the temperature distribution using equation 
(8) and, in order to flt the observed spectrum, we vary the disk inclination angle and the 
height of the wall, which are properties not constrained by other means. Thus, for each 
composition we have calculated a small grid of wall models, for Tq =130, 135, 140, 145, 150, 
155 and 160 K, cosi =0.25, 0.35, 0.45, 0.55, 0.65 and 0.75, H^^aii = 200, 220, 240, 260, 280 
and 300 (i.e., 1.77, 1.95, 2.13, 2.3, 2.48 and 2.66 AU). Assuming that gas and dust have 
the same temperature at the wall, we can evaluate a characteristic vertical gas scale height 
as hy,aii = Cs{To)/nk{Rwaii) = 1.05{To/U0Ky/^{R/10AU)^/^ , thus, the waU height we are 
trying span from 1.5 to 2.5 times hwaii- The values of H^aii/hwaii are consistent with the 



^http : //www . astro . pr inceton . edu/ draine/ dust/dustmix . html 
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height where the optical depth to the stellar radiation becomes unity times the local gas 
scale height at ~ 10 AU, in a non-truncated accretion disk. 

Table 1 summarizes the parameters that best fit the observed spectrum for each compo- 
sition and for the two different reddening laws, Draine's and Moneti's. In this table, olmgSO 
referes to Mgo.s Feo.s Si04, pyrmg 95 to Mgo.95 Feo.05 SiOa and pyrmgSO to Mgo.s Feo.2 SiOs. 
Given the uncertainties in the reddening law, dust composition and disk inclination angle, 
we can conclude that the hole has a radius R^aii between 9 and 12 AU, or Rwaii = 10.5 ± 1.5 
AU, corresponding to a maximum temperature in the atmosphere of the wall between 140 
and 150 K. If Draine's rcddcnding law is used or the silicate dust is dominated by olivine, 
the resulting wall is slightly colder and farther away than when the Moneti's law is used or 
the dust is dominated by pyroxene. In each case, a different set of disk inclination angles 
and wall heights are required to fit the observed spectrum. The inclination angle is not well 
constrained from the comparison between the synthetic spectrum, however for different com- 
positions, inclination angles, temperatures and reddening laws, we have a robust estimate of 
the height of the wall relative to its radius, H-waii/Rwaii ~ 0.22 ± 0.02. 



3.4. Limits on dust in the inner disk 

In this section, we quantify an upper limit of the dust mass in the inner hole of 
CoKu Tau/4. The dust mean temperature is calculated assuming radiative equilibrium 
between the dust grains and the stellar radiation field. We assume the inner hole, filled with 
gas and small grains, has a constant optical depth at 10 //.m (measured from the midplane 
to the surface of the disk) as a function of radius. The mass in dust in the inner disk is given 

by 

For olivine Mgo.5 Feo.5 Si04, with a amax =0.25 /im, the absorption coefficient at the 10 /xm 
peak is KiOjum = 9 cm^ g~^, assuming a silicate dust to gas mass ratio (su ~ 0.0034. Thus, 
for TiO;um=0.01, and a wall radius Rwaii ~ 10 AU, the mass in silicates in the inner hole is 
0.007 lunar masses. 

The left panel of Figure 8 shows the contribution to the SED of the dust in the hole 
(assuming a sublimation temperature Tgu — 1400 K) for different values rio^m- The temper- 
ature of the hottest dust grains inside this inner hole is such that the emissivity is larger at 
10 /im than at 20 /im, in contrast to the observed spectrum. The only way to compensate 
for this effect is to assume a particular radial dependence of TiOjum, for instance, if Tio/tm is 
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Table 1: Parameters for the wall 
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very small in the hole and increases rapidly at i? ~ -Rwaii (which is precisely how the wall is 
defined). For these small grains, we find that the upper limit in the mass of silicates in the 
inner disk is M^^/^ < 0.0007 lunar masses. A higher mass than this limit would distort the 
10 /im silicate band. 

The right panel of Figure 8 illustrates the effect of increasing the maximum grain size. 
For a fixed value of the optical depth at 10 /xm, rio^m =0.01, we show the spectrum of 
the interior for different maximum grain sizes. Since the absorption coefficient at 10 fim 
decreases with ttmax, each spectrum corresponds to a different silicate mass. If grains are 
larger, the mass limit increases. For instance, if the grains in the inner hole have amax = 10 
/i,m, then Mg°i'^ < 0.0013 lunar masses. 

There are not enough constraints on the properties of the whole disk, but speculating 
that before the gap formation it was a typical a-irradiated disk with M = 10^^ Mq yr~^ and 
a = 0.01, the mass in gas of the unperturbed inner region (for i? < 10 AU) is Mgas = 3 x 10^^ 
Mq. For the adopted sihcate dust to gas mass ratio, one would expect to find M^''^ — 250 
lunar masses in the hole. The lost mass in grains could be locked in larger bodies or could 
have been accreted by the star. In the first case, bodies would have grown to ~ 30 km to 
account for this mass, assuming that the same power law holds even for this huge size. 



3.5. The outer disk contribution 

Given the uncertainty of the fiux at 60 ^m and that there is only an upper limit for 
the fiux at 1.3 mm {Fi, < 15 mJy, Osterloh & Beckwith 1995), it is difficult to constrain the 
outer disk. 

The wall treatment assumes plane-parallel geometry. We therefore estimate that the 
heated wall region has a radial thickness AR < i^^aib where -f^waii is the thickness of the 
wall photosphere. Then, measured either radially or vertically, 

Td ~ KrfS . (21) 

If we estimate ~ 2 cm^ g~^ (see §3.3), then 

Td ~ E(cgs). (22) 

Our disk models with typical parameters for T Tauri stars (M ~ 10~^ Mgyr^^ a = 0.01) 
yield E ~ 10 g cm~^ at 10 AU, so ~ 10, consistent with our optically-thick wall treatment 
but low enough not to contribute significantly at 1.3 mm (the mass of the wall is only ~ 10"^ 
M0). Of course, the disk structure might be quite different, with a pile-up of material at the 



- 15 - 



disk wall (a ring); in this case we can limit the amount of mass present, depending upon the 
adopted opacity, assuming that the wall and the disk are optically thin, taking the Rayleigh- 
Jeans hmit, and taking the minimum temperature of the wall, ~ 85 K (equation 13) as 
the maximum temperature of the disk. 



where the adopted opacity corresponds to a maximum grain size ~ 1 mm. Smaller or bigger 
maximum grain sizes would imply a higher mass given the same flux. To better constrain the 
outer disk, it would be very useful to measure the fluxes at 70 /im and longer wavelengths. 



What is the origin of the central clearing in the CoKu Tau/4 disk? Processes such as 
grain coagulation and disk accretion are expected to play a role in preferentially reducing the 
opacity of the inner disk relative to the outer disk due to both the high surface densities and 
shorter orbital periods in this region. So, as alluded to earlier (cf. section 3.4), a possible 
explanation for the observed SED is that grain growth has led to the production of bodies 
up to 30 km in size, rendering the inner disk optically thin. Whether or not models of disk 
accretion and grain growth can reproduce in detail the observed SED of CoKu Tau/4 is an 
interesting topic for future study. 

In the context of the model presented here, we show that the observed SED of CoKu Tau/4 

is consistent with a dramatic decrement in the dust opacity at a radial distance of 10 AU. 
An interesting interpretation of this structure is that the growth of large bodies has reached 
an advanced state, producing a planet with a mass large enough to open a substantial gap, 
allow the interior material to accrete onto the star, and prevent through orbital resonances 
the progress toward the star of disk material further out. This effect has been suggested as 
the reason for central clearings in disks around T Tauri stars (e.g.. Marsh &; Mahoney 1992, 
1993) and debris disks (Jura et al. 1993; Backman, Gillett & Low 1986). We consider this 
also to be the most probable origin for the clearing in CoKu Tau/4; in fact, the case for 
planet formation is somewhat stronger than in the older systems, after considering several 
alternatives as follows. 

Radiation pressure or Poynting- Robertson effect. In an optically thin disk of 
this size, the hfetime of dust grains 10 AU from a star like CoKu Tau/4 due to these effects 
would be about 10^ years (see, e.g. Chen & Jura 2001, or Burns, Lamy & Soter 1979). In a 
disk optically thick to most of the starhght, the time scales are substantially longer (>> 1 




(23) 



4. Origin of the hole in CoKu Tau 4 



-16- 



Myr) , as the material would have to be removed an optical depth at a time. The outer disk 
of CoKu Tau/4 is quite optically thick even in the mid infrared; so are the inner disks around 
most other Classical T Tauri stars. Thus these effects are very unlikely to be important here. 
Even if they were, they would not naturally lead to a very sharp edge. 

Shcirp change in gravitational field at 10 AU (see Duschl 1989). This effect, which 
has been invoked to explain the central clearings in galactic nuclei, could be used here if the 
disk were self-gravitating at r > 10 AU; the inner part would fall in and be accreted by the 
star on the viscous time scale (about 10^ years; Quillen et al. 2004). However, the limit on 
the mass of the CoKu Tau/4 disk from millimeter-wave observations (Osterloh and Beckwith 
1995), 10^'^ Mq, is too small for the disk to be self-gravitating. 

Radial grain segregation. In debris disks around A stars young enough for some gas 
still to be present, it has been predicted that the combined effects of radiation pressure and 
gas drag would lead to central clearings in the distribution of small grains (Takeuchi and 
Artymowicz 2001). The radiative requirements are similar to those of radiation pressure and 
Poynting-Robertson drag, though; in the present case the disk is too optically thick, and the 
star too cool, for this process to produce a central clearing within 1 Myr. 

Stellar companion. Were CoKu Tau/4 to have a stellar or planetary companion, 
a clear gap would be produced on timescales shorter than the 1 Myr stellar age, if the 
accretion timescale were sufficiently short, and if the companion were sufficiently massive 
both to open a gap in the disk (i.e., if its Hill radius were to exceed the original half- 
density thickness of the disk), and to prevent migration. Accretion of the inner-disk material 
would take place on the viscous timescale of 10^ years, which is indeed sufficiently short. A 
terrestrial-size planet would be unable to produce a gap or prevent migration. A stellar 
companion could do both. No such companion has yet been observed. The visible and 
infrared properties of CoKu Tau/4 are accurately consistent with a single Ml. 5 pre-main- 
sequence star; a companion star is constrained to have luminosity of at most a small fraction 
of that of the Ml. 5 star, for which the luminosity is already small on the scale of stars. Thus 
it is not impossible that the central clearing has been produced and maintained by a very 
low-mass star which has so far escaped detection. The range of stellar masses for this to work 
is so narrow, however, that this seems less likely an option than a substellar companion. 

A companion in a wide range of masses in the giant planet or brown dwarf regime 
would produce such an empty, sharp-edged clearing as we see. Upon formation of the 
companion, material further out than the strongest orbital resonances would be held up by 
frequent collisions at these resonances until the mass held up were similar to the mass of the 
companion. Material inside these resonances would simply be accreted on the viscous time 
scale. Quillen et al. (2004) have recently simulated the dynamics of a CoKu Tau/4-like disk 
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that undergoes giant-planet formation, and found that clearings like the one observed were 
produced for companions as small as 0.1 Jupiter masses. 

Observations can be used to test the scenario proposed in this paper to interpret the 
spectrum of CoKu Tau/4. For instance, molecular hne observations could probe whether the 
outer disk is gaseous; continuum millimeter observations can be used to quantify the outer 
disk mass; spatially resolved images could confirm that there is a sharp decline in the dust 
column density around 10 AU and might reveal whether the disk has the non-axisymmetric 
structure that is expected to be induced by the planet (e.g., Quillen et al. 2004). 

5. Discussion and Sumniciry 

The spectrum of CoKu Tau/4 is consistent with that expected from a disk truncated 
at an inner radius ~ 10 AU. No detectable material is in the inner hole; the spectrum 
indicates at most 0.0007 lunar masses of small silicate grains in the inner hole region. Even 
though CoKu Tau/4 is much younger than TW Hya, the hole in its inner disk is larger and 
apparently more evacuated of dust and gas (judged by the lack of near-infrared dust emission 
and gaseous accretion), and its outer disk appears to be much less massive. This comparison 
clearly indicates that age is not the sole parameter in determining disk evolution. That other 
parameters influence the properties of pre-main sequence disks has been clear for some time. 
For example, nearly half the stars in Taurus lack inner disks, some (but not necessarily all) 
of which can be explained by the presence of perturbing binary companions (e.g., Mathieu 
et al. 2000, and references therein). 

Because it seems likely that even binary systems once possessed accretion disks (see, 
e.g.. Bate, Bonnell, & Bromm 2003), there must be a wide range of evolutionary timescales 
for protoplanetary disks. The reasons for this range of timescales are not clear, but we 
conjecture that initial angular momentum could play a crucial role; disks of similar mass 
but much smaller initial radii are more likely to fragment and to coagulate faster (given 
higher initial disk surface densities). This conjecture predicts that the outer radius of the 
CoKu Tau/4 disk is small; observational efforts to explore this possibility should be made. 
In any case, given the low disk mass estimate, CoKu Tau/4 may be a system observed just 
barely before the initial disk material disappears. 

Finally, to summarize our results, we find that an inner disk wall at i? ~ 10 AU with a 
half-height i^waii ~ 2 AU , is responsible for the observed excess between 8 and 30 /xm. This 
inner wall has an optically thin atmosphere, closer to the star, with a temperature decreasing 
with radius. The highest temperature of the wall is ~ 145 K, but decreases rapidly to ~ 85 
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K in the deepest layers. Both the 10 and 20 //m bands emerge from the optically thin hotter 
layers. The spectrum for A > 20 urn emerges from the optically thick layers of the wall. At 
even longer wavelengths the whole wall becomes optically thin, and it is expected that the 
outer disk dominates the SED. Using an upper flux limit at 1.3 mm (Osterloh & Beckwith 
1995), and assuming that the wall is so optically thin at this wavelength that it does not 
contribute to the flux, we estimate the mass of the outer disk to be less than 0.001 for 
a maximum grain size of 1 mm. The central star dominates the spectrum at A < 8 /xm. 

From an analysis of the 10 fim feature, wc find that most of the dust in the wall should 
be glassy Fe Mg olivine and/or pyroxene, with no evidence of a crystalline component. We 
believe this represents the initial, unprocessed state of the dust in the disk that formed 
around CoKu Tau/4 . 



6. Appendix 

The coordinate system x, y and z is centered on the star, the z axis is the disk rotation 
axis, and the plane {x, y) is the disk midplanc. Wc use also cylindrical coordinates R, 9 and 
z. The wall is a cylinder with radius -Rwaii and total height 2if„aii- Points on the wall surface 
have coordinates z,x — RaCOsO and y = RasinO. We use another coordinate system X, 
Y, Z, also centered on the star, where Z is along the fine of sight and {X, Y) is the plane 
of the sky. Both systems coincide when the disk is pole-on. The transformation between 
coordinate systems is 

x^X (24) 
y = Y cos i — Zsini (25) 
z = Y smi + Z cos i (26) 

or 

X = x (27) 
Y — y cos i + zsini (28) 
Z — z cos i — y sin i (29) 

The normal to the cylinder farthest wall (which is the one the observer will see) is 
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n — —R — — cos Ox — sin Oij (30) 
thus, the cosine between the hne of sight 

Z = — sini ^ + cosi i, (31) 
and the normal to the visible surface of the wall is 

cos Q — Z ■ h — sini sin ^ (32) 

A surface element of the visible area is dA — dX dY — —Rd sin 9 dO dY . To find the 
limits of the visible surface, we have to define two ellipses (given by the projections of the 
upper and lower tops of the wall; see Figure 9). For the upper eUipse 

^ = ± cos i a/i - X Vi?2 + ^ sin i (33) 

Rd ^ Rd 

For the lower ellipse. 



± cos i Jl-X'^/Rl - — sin i (34) 



Rd ^ Rd 

Which can be written also in terms of 9, using X — x — RaCOsO, 

= ± cos i sin 6* + — sin i (35) 

Rd Rd 
^down I . . ^ H , . fnn\ 

— — = ± COS ismO — — sm i (36) 

Rd Rd 

Both ellipses intersect at critical angles: 6c and ir — 9c, where 9c is given by 

H H 
— cos i sin 9c + — sin i — cos i sin 9c — — sin i (37) 

Rd Rd 

sin 9c = — tan i — 5 (38) 

Rd 



Depending on the inclination angle i, both ellipses can intersect if 5 < 1 or not if 5 > 1 
(see Figure 9). U S > 1, then the integration region is defined by the upper ellipse and 
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Y~p<Y< Y^, < 9 < TT. If5<l, there are two regions: if < 6* < sin ^ 5 then 
yup<Y< but if sin-i S<9<n/2 then Y-^<Y< 
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Fig. 1. — Spectrum of Co Ku Tau 4 corrected using reddening laws by Draine (2003) and 
Moneti et al. (2001) assuming Ay=3. The IRS fluxes are shown by error bars wich represent 
the observed error, estimated as half of the absolute value of the difference in flux from the 
spectra independently obtained at the two nod positions of the telescope for each order of 
Short-Low and Long-Low. 
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Fig. 2. — Spectral energy distribution of CoKu Tau/4 (squares from Strom et al. 1989; 
circles from 2MASS; IRS portion from Forrest et al. 2004) corrected using the reddening 
law by Moneti et al. (2001) with Ay — 3, compared to the median SED of Classical T Tauri 
stars in Taurus (triangles and dashed line). The error bars of the median points represent 
the quartilcs (D'Alessio et al. 1999). We also show a synthetic spectrum for the central star 
taken from Bruzual & Chariot (1993), for % = 3720 K and = 1.9 R©. 
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Fig. 4. — Radial distribution of temperature for wall models whose synthetic spectrum fits 
the observed CoKu Tau/4 spectrum. Both models have Tq =145 K, one contains Mgo.5 Feo.5 
olivine (solid line) and the other one has Mgo.s Feo.2 pyroxene (dashed line), to illustrate the 
effect of the different silicate compositions on the temperature distribution. 
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Fig. 5. — Spectra of wall models calculated using diferent optical constants: silicates from 
Draine(2003) (dotted line) for H^aii =280 R^,; glassy silicate of mean cosmic composition 
M go. 50 Feo.43 Cao.o3 AI0.04 SiOs (dashed line) with optical constants from Jager et al. (1994) 
and H^^aii =240 i?* ; glassy bronzite (solid line) with optical constants from Dorschner et al. 
(1988) and Hu,aii = 240 i?*; oxygen rich circumstellar sihcates (dot-dashed hne) with optical 
constants from Ossenkopf et al. (1992) and H^aii — 320 R^. In all these models Tq =145 
K and cosi = 0.55. The synthetic spectra are compared to the spectrum of CoKu Tau/4 
corrected using Moneti's reddening law. 
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Fig. 6. — Examples of wall model spectra (solid lines) that fit to the observed spectrum of 
CoKu Tau/4 corrected by the reddening laws of Draine (2003) (model parameters: Tq = 
140 K, cosi =0.45 and Hu,aii =280 R^) and Moneti et al. (2001) (model parameters: Tq = 
145 K, cosi =0.45 and Hy^aii =260 R^). The dust consists of small grains {amax — 0.25 //m) 
of glassy ohvine, with 50 % Fe and 50 % Mg, with a small amount of organic and troilite 
grains (see §3.3). 
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Fig. 7. — Example of wall models spectra (solid lines) that fit to the observed spectrum of 
CoKu Tau/4 corrected by the reddening laws of Draine (2003) (model parameters: Tq = 
140 K, cosi =0.45 and Hu,aii =240 R^) and Moneti et al. (2001) (model parameters: Tq = 
145 K, cosi =0.45 and Hy^aii =240 R^) . The dust consists of small grains {amax — 0.25 //m) 
of glassy pyroxene, with 20 % Fe and 80 % Mg, with a small amount of organic and troilite 
grains (see §3.3). 
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Fig. 8. — Spectrum of dust in the inner disk compared to the IRS spectrum of Co Ku Tau 
4 corrected using the reddening law by Moneti et al. 2001. The adopted dust is composed 
by glassy olivine, with 50 % Fe and 50 % Mg. Left panel: for a^ax = 0.25 jim. and each 
curve corresponds to a different total sihcate mass M^^'.f^ = 0.00015, 0.0004, 0.0007, 0.0015, 
0.004, 0.007, 0.015 and 0.04 lunar masses (corresponding to tio^^ =0.0002, 0.0005, 0.001, 
0.002, 0.005, 0.01, 0.02, 0.05, from bottom to top). Right panel: for tiq^^ = 0.01; each 
curve corresponds to a different maximum grain size amax — 0.25 (solid line), 10 (long- 
dashed hue) (corresponding to M^"^^ — 0.007 and 0.013 lunar masses, respectively). Biger 
maximum grain sizes, for the same optical depth at 10 ;um, produce a similar spectrum than 
the case for a.max = 10 ;U m, but for a different total dust mass. For instance, a^ax = 1 cm 
corresponds to 0.36 lunar masses, a^aax = 1 ni, to 3.6 lunar masses. In each of these cases, 
dmin = 0.005 iim and the exponent of the size distribution is p = 3.5. 
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Fig. 9. — Schematic representation of the visible surface of the wall (dashed surface) as seen 
by the observer for two inclination angles, one corresponding to 5 < 1 (the star is visible) 
and the other one to 5 > 1 (the star is invisible). 



